Universality in a Class of Q-Ball Solutions: An Analytic Approach 
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The properties of Q-balls in the general case of a sixth order potential have been studied using 
analytic methods. In particular, for a given potential, the initial field value that leads to the soliton solution 
has been derived and the corresponding energy and charge have been explicitly evaluated. The proposed 
scheme is found to work reasonably well for all allowed values of the model parameters. 

1 Introduction 

A scalar field theory with a spontaneously broken U{1) symmetry may contain stable nontopological solitons, 
the so-called Q-balls [l]-[2]. Q-balls are coherent states of complex scalar fields that carry a global U{\) charge 
and can be understood as bound states of scalar particles which appear as stable classical solutions carrying 
a rotating time dependent internal phase. They are characterized by a conserved nontopological charge Q 
(Noether charge) that ensures existence and stability [3]- [4]. 

The concepts associated with these solutions are quite general and occur in a wide variety of physical 
contexts [5] . Q-balls are allowed in supersymmetric extensions of the standard model that allow flat directions 
in the scalar potential. Flat directions in the Minimal Supersymmetric Standard Model (MSSM) [6] have 
been shown to exist. The conserved charge is associated with the U{\) symmetries of baryons and leptons, 
while, the relevant U{1) fields correspond to either squark or slepton particles. Thus, Q-balls can be thought 
of as condensates of a large number of either squark or slepton particles which can affect baryogenesis via 
the Affleck-Dine mechanism [7] during the post-inflationary period of the early universe. The Q-ball stability 
is cosmologically important since if stable Q-balls are formed in the early universe they can contribute to 
its dark matter content. These can be huge balls with charges of order 10^°; however, very small Q-balls 
can also be considered as dark matter constituents [8] . Decaying Q-balls can also be of crucial cosmological 
significance. If Q-balls decay after the electroweak phase transition, they can protect baryons from the 
erasure of baryon number due to sphaleron transitions. Furthermore, the Q-ball decay may contribute to 
dark matter production in the form of the lightest supersymmetric particle, explaining the baryon to dark 
matter ratio of the universe [9] . 

Consider the U{1) Goldstone model Lagrangian describing a single complex scalar field </) in three 
spatial dimensions given by 



?7(|(/)|) has a single minimum at = which is equivalent of stating that there is a sector of scalar particles 
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c = -d^cpd^^-um). 



(1) 



(mesons) carrying U{1) charge with mass equal to J ^U"{0). The corresponding energy functional is given 




while the conserved Noether current (due to the global U{1) symmetry) is 




(3) 



with charge given by 




(4) 



The stationary Q-ball solution can be obtained by assuming that 



- e-*/(r) (5) 
where /(r) is a real radial profile function that satisfies the ordinary differential equation 

+ + (6) 

dr^ r clr 

with boundary conditions /(oo) = and f'{0) = 0. In each case, the effective potential is defined as 
Ueffif) = iJ^ p 12 — U{f), while, the existence of Q-ball solutions leads to constraints on the potential U{f) 
and the frequency w: (i) The effective mass of / must be negative, so, by assuming that C/(0) = U'{0) ~ 
and U"{0) = u)"^ > one can deduce that iv < U!+. (ii) The minimum of U{f)/ p must be attained at 
some positive value of / (say < /o < oo) and existence of the solution requires that lo > U- where 
uP_ = 2t/(/o)//o • Hence, Q-balls exist for all u) in the range uj^ < < cj+. 

The charge (4) and energy (2) of a stationary Q-ball solution (5) take the simple form 



Jr^fir)dr (7) 
E = Ito^Q^ +4n [ (^ + U{f)]Pdr. (8) 



Q = Attlj 



2 J \ 2 

Numerical and analytical methods have shown that when the internal frequency is close to the minimal value 

cij_ the profile function is almost constant, implying that the charge (7) is large (thin- wall approximation). 
On the other hand, when the internal frequency approaches the maximal value uj-\. the profile function falls 
off very quickly (thick- wall approximation). In the thick- wall approximation the behavior of Q depends on 
the particular form of the potential and the number of dimensions [11]. 

The choice of the potential is not unique since the only requirement is that the ratio U{f)/ p has a 
local minimum at some value of / different from zero. There are several natural types to be considered, two 
of which are shown below: 

I: U{f) = ^f-X,f + ^,f, 

II: [/(/) = msf{l-K\og(Xs.f))+fisf''. (9) 

In each case, two of the parameters can be removed by rescaling, thus, the potentials of type I have one 
free parameter while potentials of type II have two for fixed p. The type I potential is the simplest allowed 
one which is a polynomial in p, while type II mimics the D-flat direction in MSSM. Here, p > 6 is some 
integer that ensures the growth of the potential for large /, but does not destroy the flatness property for 
intermediate values of /. None of these types are the kind which might be associated with a renormalizable 
quantum field theory, but, they arc typical of effective theories incorporating radiative or finite temperature 
corrections to a bare potential. In this paper, we shall be dealing with type I potentials only, for which 

a;+ = mi and a;_ = \f^^^^^^ ■ Stable Q-balls exist for y^ ^™ 2m\~ ~ < ijJ < mi. The equation of motion 
(6) is 

^^^ll = „^f-,,,f^,,,f. (10) 

where we have defined = m\ — uP' . 

Q-balls can be studied either analytically [1. 11, 14], or numerically [3]- [13]. In a recent work [15, 16], 
we employed analytic arguments in order to construct an approximate profile function of the symmetrized 
Woods-Saxon type, expected to be valid in the thin-wall limit. This approximate profile led to an explicit 
energy-charge relation which, to our surprise, was found to yield valid results for a region far exceeding the 
expected limits of the thin- wall approximation. The calculation was carried out for a specific form of type 
I potential i.e. Ai 2, /xi = 1 and m = 2. In this paper, the aforementioned work is being extended to 
include the general case of type I potentials (where the parameters are taken to be arbitrary) in hoth regimes, 
thin- wall and thick- wall respectively. This way, the soliton energy and charge as well as the initial field value 
/(O) can be accurately derived as functions of a for the whole allowed range of the model parameters. 



2 Potential Energy 



From a mechanical point of view, a Q-ball solution describes the motion of a particle moving with friction 
in the potential 

C^e//(/) = -^aV' + Alr-Ml/^ Ai,Mi>0. (11) 

The corresponding equation of motion is given by (10). Upon rescaling r — > r/a and letting /(r) = f{0)ip{ar) 
for ip{0) = 1, equation (10) transforms to 

(filj 2dilj ^ 
ar^ r dr 

where A = ^^^Ji^ and jj, = and the effective potential (11) becomes 

C^W = -^^' + AV'-/^^'- (13) 
After some algebra, it can be shown that equation (12) implies that 

r'-""'(^)'*=i^r'-"""(-5*'^ 

For the special cases n = 3 and n = 0, we get that 



llr^if)'' = 3/^^(-i^^ + A^^-.^«).. (15) 



2j\-'(^^ydr = -iv'(0)2 + AV(0)^-M^(0)'^ 

= C/(l). (16) 

The first condition is Derrick's (or virial) theorem which states that, for a three dimensional model, the 
kinetic energy equals three times the potential energy; the second condition describes the energy dissipation 
of the mechanical system. 

The mechanical analogue imposes several constraints on U{1). Consider a particle initially located 
at the point ip{0) — 1 which starts rolling down the potential wall to eventually stop at the point 1/1(00) = 0. 
(Here r corresponds to the time variable). Then, the potential energy (13) at the origin i7(l) = — 5+A — /U 
has to be consumed by the friction term therefore U{1) has to be bounded both from below and above: 
since for U{1) large the particle will overshoot the top point while for U{1) small it will not reach it. 

The initial potential energy U{1) has to be positive implying that 

M < A - ^ (17) 

and, also, U'{1) has to be positive (attractive "force") leading to 

^ < i (4A - 1) . (18) 

These conditions arc satisfied for A > 1 when (18) holds while for 1/2 < A < 1 when (17) holds. Note that 
for A < 1 no Q-ball solutions can be found. Depending on the shape of the potential (i.e. the actual values 
of A and /i) the following two distinct cases occur: 

/. Thin-Wall Approxim,ation: In this case, V'(O) lies near the maximum of the effective potential 
which is deep. Then U'{1) = — 1 + 4A — 6/i « has to be positive and close to zero (slow roll) while U"{1) 
is negative (convex region) and of order unity. The Q-ball solutions lie approximately on the line 

M-^(4A-1) (19) 



and the initial potential energy depends linearly on A: 



C/(l) = i(A-l), A>1. (20) 

//. Thick-Wall Approximation: Here, the potential is shallow and the maxima arc high up while 
U'{1) is positive and large (concave region) and fx is small. The initial potential energy U{1) increases with 
A and reaches its maximal value when — > 0. This value can be determined numerically and is found to 
be: Xmax — 4.70137. The actual functional dependence of U{1) on A however, is unknown and needs to be 
determined. 

To that end, let us assume a leading power law behaviour 

C/(l) kA". (21) 
Then, (13) implies that the Q-ball solutions lie on the line 

/i = -^ + A-KA" (22) 
which determines k in terms of Xmax (since n{Xmax) = 0): 

^^2W^_ (23) 

^max 

In order to determine the value of n in (21), we require that the transition from the thick-wall region 
to the thin-wall one be smooth implying that, U'{1) is small for small values of A (~ 1), i.e. 

C/'(1)=2(1-A)+6kA"«0. (24) 
1 

Its minimum value occurs at the point Ajnin = (3nfi;)i"". Substituting Amin in (24) one gets that U'{1) 
vanishes for n = 2.82187 in which case Amin = 1.54889 (in agreement to our hypothesis). Assuming that 
n can only take integer or half integer values (otherwise the structure of the model on the complex plane 
would be very complicated). Figure [1] shows that for n = 3, the transition is indeed satisfactorily smooth. 
For this case, we get that k = ■ This way, we get the /x(A) relation for the thin-wall (19) and thick-wall 
(22) limits depicted in Figure [2] against numerical data. 



3 Properties of Q-balls 

Using the rescaling formulae of the previous section, the a dependence of the functions E, Q and /(O) can 
now be analytically determined. The initial field value /(O) that leads to a Q-ball solution is given by 

Figure [3] depicts predicted values of /(O) in terms of a against numerically calculated ones for sample values 
Ai = 2, /xi = 1, m = 2 (as in [16]). 

The corresponding charge (7) and energy (8) functionals are given by: 

^ ^ 4m. A r^2^2^^ (26) 
a Ai Jq 

E = ^±^Q + 47ra^ r (l^'^-XiP^+fiijAr^dr- (27) 
2uj Xi Jo \2 ) 



It was shown in [16] that a trial function satisfying the boundary conditions and having the right 
asymptotic behavior, is the symmetrized Woods-Saxon profile 

= ' (28) 

^^1 + ci cosh (or) 




Figure 1: The potential energy gradient in terms of A for different values of n. 




Figure 2: The /i(A) relation for n = 3 against numerical data. 




Figure 3: Predicted initial field values /(O) against numerical data. 

This function satisfies the exact differential equation 

^"-^ V + ^^(l-c?)0^ (29) 



and the approximate one 



Here, c, ci and h are arbitrary parameters which need to be determined in order that (28) fits the exact profile 
function ^{r) in the best possible way. Note that only two of the three parameters in (28) are independent 
since the initial condition 0(0) = 1 implies that c = ^/l + ci. 
First note that (28) satisfies the following relations: 

/)2 1 /)2 1 1,2 

a4 2/2,2 ^ J.2 

aci 

4>^ = c'4>'' + 2c'c^^<jP + \d'cl^cf' (31) 
ac\ 2 aci 

in terms of which the charge (26) and energy (27) functional can be explicitly evaluated: 
^ 4tvlj a . 

Q = —^^^0 (32) 



a Ai 363 




dci 

(33) 



Here io is a function of ci given by 



io (ci) 



arccosh ( — | tt^ + arccosh ( — ) 

\ClJ \ClJ 



for ci < 1. (34) 



Note that Ai, fii, m are fixed external parameters, A varies continuously between 1 and 4.70137 and /x is a 
known function of A given by (19) or (22) depending on the regime. 

The profile parameters b and ci can be determined by imposing the conditions (15) and (16) on 
4>{r), that is 



1 /"^ 



dr sinh^ r 



^ = Uil) (35) 

r (1 + ciCOshr)^ 

^_...,^._!M),,(.^._5M_v)c*4(.c'4(:-.n)f = 0. (3e, 



Equations (35) can be written in a compact form 




(37) 



where 



in terms of which (36) becomes 



dr sinh^ r 



ii= [ — '^^^^ ^ (38) 

io r (1 + ci coshr)'* 



1,2 4 1 U(l)\ . / 4 1 U(l) , n\ dio 41 , (l-cl) U(l)\ d'^io „ 
_ _ Ac^ + _ _ J_i j + ^2,c^ _ _ ^ _ Ac^ j ^ + ^ (^,0^ + 

(39) 

Equation (39) determines ci in terms of A and and is depicted in Figure [4]. Note that the values of ci and 
b obtained this way are universal since they do not depend on the geometrical parameters of the potential. 
Nevertheless, the values of the energy and charge do depend on the specific form of the potential i.e. the 
values of Ai, /xi and m. These values tend to infinity in both limits since (i) in the thin- wall limit Ci ^ 
while (ii) in the the thick-wall limit b 0. 

Figures [5] and [6] presents the a dependence of the charge (32) and energy (33) against values 
obtained numerically, for Ai = 2, /xi = 1, m = 2. It is interesting to realize that the range of validity of the 
thick-wall approximation is very wide and gives satisfactory results even in the thin-wall region. It appears 
that in this class of theories the Q-balls prefer to be "thick" . 



4 Conclusions 

In this paper we have extended our earlier work [16] and investigated phenomenologically relevant properties 
of Q-balls in a universal way. In particular, we have addressed the following problem: Given the geometrical 
characteristics of the scalar potential find the initial field value that will lead to a Q-ball solution as well 
as the corresponding charge and energy of the soliton. For a particular class of potentials (sixth order 
polynomials), after parameter rescaling, the problem can be tackled in a universal way. The Q-ball profile 
can be accurately approximated by means of a two-parameter symmetrized Woods-Saxon function which 
can be analytically calculated in all cases. This scheme is found to yield satisfactory results in the whole 
parameter region so that we do not have to rely on approximations like thin-wall or thick-wall. The soliton 
energy and charge can subsequently be analytically calculated and compared against numerically calculated 
values taken from an earlier work [15]. 

We believe that a similar line of argument can be applied to study the profile function and the 
energy-charge dependence in all types of potentials. 
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Figure 4: The Ci(A) relation obtained from (39) 




Figure 5: The a dependence of Q given by (32) 




Figure 6: The a dependence of E given by (33). 
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